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Quantum networks are composed of nodes which can send and receive quantum states by exchang- 
ing photons Their goal is to facilitate quantum communication between any nodes, something 
which can be used to send secret messages in a secure way @,|3|], an d to communicate more efficiently 
than in classical networks These goals can be achieved, for instance, via teleportation [jj. Here 
we show that the design of efficient quantum communication protocols in quantum networks involves 
intriguing quantum phenomena, depending both on the way the nodes are displayed, and the en- 
tanglement between them. These phenomena can be employed to design protocols which overcome 
the exponential decrease of signals with the number of nodes. We relate the problem of establishing 
maximally entangled states between nodes to classical percolation in statistical mechanics [6|, and 
demonstrate that quantum phase transitions [7J can be used to optimize the operation of quantum 
networks. 
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Introduction. The future of quantum communication 
will be based on quantum networks (cf. [l|, [H, [^, Ell El, 
[HI EH, EH ) > where different nodes are entangled, leading 
to quantum correlations which can be exploited by per- 
forming local measurements in each node. For instance, 
a set of quantum repeaters 10] can be considered as a 
simple quantum network where the goal is to establish 
quantum communication over long distances. In order to 
optimize the operation of such a network, it is required 
to establish efficient protocols of measurements in such 
a way that the probability of success in obtaining max- 
imally entangled states between different nodes is maxi- 
mized. This probability may behave very differently as a 
function of the number of nodes if we use different proto- 
cols: in some cases it may decay exponentially, something 
which makes the repeaters useless, whereas for some pro- 
tocols it may decay only polynomially, something which 
would make them very efficient. 

A general network may be characterized by a quantum 
state, p, shared by the different nodes. The goal is then: 
given two nodes, A and B, find the measurements to be 
performed in the nodes, assisted with classical commu- 
nication, such that A and B share a maximally entan- 
gled state, or singlet, with maximal probability. We call 
this probability the singlet conversion probability (SCP). 
This, or other related quantities like the localizable en- 
tanglement [HI, [l6[ , can be used as a figure of merit to 
characterize the state p and therefore the performance of 
the quantum network. Here, we focus on the SCP be- 
cause of its operational meaning. These quantities can- 
not be determined in general, given that they require the 
optimization over all possible measurements in the dif- 
ferent nodes, which is a formidable task even for small 
networks. 

In this work we concentrate in some particular quan- 
tum networks which, despite its apparent simplicity, con- 



tain a very rich and intriguing behavior. The simplifica- 
tion comes from two facts (see Fig. 1): first, the nodes are 
spatially distributed in a regular way according to some 
geometry. Second, each pair of nodes are connected by a 
pure state \ip) <E <C d <£> <C d . Up to local change of bases, 
any of these states can be written as [TtJ 
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where A^ are the (real) Schmidt coefficients such that 
Ai > A2 > . . . > Xd > 0. This configuration reminds 
of the states underlying the so-called projected pair en- 
tangled states [lH , and thus we call these networks pair- 
entangled pure networks (PEPN). For these geometries, 
we first introduce a series of protocols which are closely 
related to classical percolation [|| , a concept that appears 
in statistical mechanics. We then determine the optimal 
protocols for several 1-dimensional (ID) configurations, 
where some counterintuitive phenomena occur. We use 
these phenomena to introduce various protocols in more 
complex 2-dimensional (2D) configurations. 

We show that these new protocols provide a dramatic 
improvement over those based on classical percolation, 
in the sense that one can obtain perfect quantum com- 
munication even though the percolation protocols give 
rise to an exponential decay of the success probability 
with the number of nodes. In fact, we will argue that 
there exists a quantum phase transition in the quantum 
networks which may be exploited to obtain very efficient 
protocols. Thus, this work opens a new set of problems 
in quantum information theory which are related to sta- 
tistical physics, but pose completely new challenges in 
those fields. As opposed to most of the recent work on 
entanglement theory, which has been devoted to using 
some of the tools developed so far in quantum informa- 
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tion theory to analyze problems in statistical mechanics 
[H> [3 HU IH|, the present work makes a step in the 
converse direction. 

Classical entanglement percolation. A first natural 
measurement consists of all the pairs of nodes nodes lo- 
cally transforming their states into singlets with optimal 
probability, p ok . Recall that the SCP for a state (p} is 
known to be equal to p ok = min(l, 2(1- Ai)) [13]. Then, 
a perfect quantum channel between the nodes is estab- 
lished with probability p ok , otherwise no entanglement 
is left. This problem is equivalent to a standard bond 
percolation situation [H, where one distributes connec- 
tions among the nodes of a lattice in a probabilistic way: 
with probability p an edge connecting a pair of nodes is 
established, otherwise the nodes are kept unconnected. 
We call this measurement strategy classical entanglement 
percolation (CEP). In bond percolation, for each lattice 
geometry there exists a percolation threshold probabil- 
ity, pth, such that an infinite connected cluster can be 
established if and only if p > p t h (see also Table 1). The 
probability 9{p) that a given node belongs to an infinite 
cluster, or percolation probability, is strictly positive for 
P > Pth, and zero otherwise (in the limit of an infinite 
number of nodes). Then, the probability that two given 
distant nodes can be connected by a path is distance inde- 
pendent and given, correspondingly, by 9 2 (p) for p > p t h] 
for p < pth this probability decays exponentially with the 
number of nodes, N. 

The threshold probabilities define a minimal amount 
of entanglement for the initial state such that CEP is 
possible. In the case of ID chains, see Fig. [2] percolation 
is possible if, and only if, p = 1. Therefore, the SCP 
decays exponentially with TV unless the states are more 
entangled than the singlet, in the sense that p ok = 1. In 
a square 2D lattice, the entanglement threshold derived 
from percolation arguments, see Table HI is p ok = 2(1 — 
Ai) = 1/2. 

It is natural to wonder whether CEP is optimal for 
any geometry and number of nodes. And if not, to see if, 
at least, it predicts the correct decay of entanglement in 
the asymptotic limit. Next, we show that, for ID chains, 
although CEP is not optimal for some finite N, it gives 
the right asymptotic behavior. Moving to 2D networks, 
we prove that CEP is not optimal even in the asymptotic 
case. Thus, the problem of entanglement distribution 
through quantum networks defines a new type of critical 
phenomenon, that we call entanglement percolation. 
ID Chains The scenario of a ID chain configuration, see 
Fig. [2] consists of two end nodes connected by several 
repeaters [Hj . As said, all the bonds are equal to \ip). 

We start out with the case of qubits, d — 2. A surpris- 
ing result already appears in the first non-trivial situa- 
tion consisting of one repeater. An upper bound to the 
SCP in the one-repeater scenario is obtained by putting 
nodes A and R\ together, which implies that the SCP 
cannot be larger than p ok . This bound can indeed be 
achieved by means of a rather simple protocol involving 
entanglement swapping 25] at the repeater. However, 



if CEP is applied, the obtained SCP is simply (p ok ) 2 . 
This proves that CEP is not optimal already for the one- 
repeater configuration. We find quite counter-intuitive 
that the intermediate repeater does not decreases the op- 
timal SCP. This behavior, however, does not survive in 
the as ymp totic limit. In this limit, the so-called concur- 
rence [241 ] . another measure of entanglement, decreases 
exponentially with the number of repeaters, unless the 
connecting states are maximally entangled (see Meth- 
ods). The exponential decay of the SCP automatically 
follows. 

Most of these results can be generalized to higher di- 
mensional systems, d > 2. For the one-repeater configu- 
ration, the SCP is again equal to p ok . It suffices to map 
the initial state into a two-qubit state, without changing 
the SCP, and then apply the previous protocol. Mov- 
ing to the asymptotic limit, an exponential decay of the 
SCP with N can be proven in the scenario where the 
measurement strategies only involve one-way communi- 
cation: first, a measurement is performed at the first 
repeater. The result is communicated to the second re- 
peater, where a second measurement is applied. The re- 
sults of the two measurements are communicated to the 
third, and so on until the last repeater, where the final 
measurement depends on all the previous results. 

Putting all these results together, a unified picture 
emerges for the distribution of entanglement in ID 
chains: despite some remarkable effects for finite N, 
the SCP decreases exponentially with the number of re- 
peaters whenever the connecting bonds have less entan- 
glement than a singlet. The CEP strategy fails for some 
finite configuration, but predicts the correct behavior in 
the asymptotic limit. 

2D Lattices. The situation becomes much richer for 
2D geometries. First, we consider finite 2D lattices. The 
non-optimality of CEP can be shown already for the sim- 
plest 2x2 square lattice and qubits. Consider the two 
non-neighboring sites in the main diagonal of the square. 
The SCP obtained by CEP is 1 - (1 - (p ok ) 2 ) 2 . By con- 
catenating the optimal measurement strategy for the one- 
repeater configuration, the SCP is 1 — (1 — p ok ) 2 . How- 
ever none of these strategies exploits the richness of the 
2D configuration. Indeed, one can design strategies such 
that a singlet can be established with probability one 
whenever \ip) satisfies 1/2 < Ai J$ 0.6498. Thus, there 
are 2D network geometries where, although the connec- 
tions are not maximally entangled, the entanglement is 
still sufficient to establish a perfect quantum channel. 

Let us now see whether the thresholds defined by stan- 
dard percolation theory are optimal for asymptotically 
large networks. In the next lines, we construct an ex- 
ample that goes beyond the classical percolation pic- 
ture, proving that the CEP strategy is not optimal. The 
key ingredient for this construction is the measurement 
derived above for the one-repeater configuration, which 
gave raise to a SCP equal to p ok . Our example consid- 
ers a honeycomb lattice where each node is connected by 
two copies of the same two-qubit state \<p) = l^)® 2 , see 
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Fig. [3la. If, as above, the Schmidt coefficients of the 
two-qubit state are Xi > A 2 , the SCP of \<p) is given by 
p ok = 2(1 — Xf). We choose this conversion probability 
smaller than the percolation threshold for the honeycomb 
lattice, which gives 

Al = ^ +Sin (S)"°- 82 - (2) 

So, CEP is useless. Now, half of the nodes perform the 
optimal strategy for the one-repeater configuration, map- 
ping the honeycomb lattice into a triangular lattice, as 
shown in Fig. [3j The SCP for the new bonds is exactly 
the same as for the state l^a), that is 2X 2 . This prob- 
ability is larger than the percolation threshold for the 
triangular lattice, since 

2A 2 = 2 (l - * °' 358 > 2sin (£) • 

(3) 

The nodes can now apply CEP to the new lattice and 
succeed. Thus, this strategy, which combines entangle- 
ment swapping and CEP, allows to establish a perfect 
quantum channel in a network where CEP fails. 
Conclusions. We have shown that the distribution 
of entanglement through quantum networks defines a 
framework where statistical methods and concepts, such 
as classical percolation theory and beyond, naturally ap- 
ply. It leads to a novel type of quantum phase transi- 
tions, that we call entanglement percolation, where the 
critical parameter is the minimal amount of entanglement 
necessary to establish a perfect distant quantum channel 
with significant (non-exponentially decaying) probability. 
Further understanding of optimal entanglement percola- 
tion strategies is necessary for the future development 
and prosperity of quantum networks. 

METHODS 
ID chains 

We start by showing that the concurrence decays ex- 
ponentially with the number of nodes in a ID chain of 
qubits when the connecting states are not maximally en- 
tangled. Recall that, given a two-qubit pure state \ip) = 
Ylij > its concurrence reads C(<p) = 2|det(T)|, 
where T is the 2x2 matrix such that (T)ij = tij. 

When considering the repeater configuration, the max- 
imization of the averaged concurrence turns out to be 
equal to 

C N =supy)2\det(cp 1 M ri <p2...M rN <p} T+1 )\. (4) 

M — 

r 

Here M. briefly denotes the choice of measurements, 
while ipk represent the 2x2 diagonal matrices given by 
the Schmidt coefficients of the states \(pk)- M rk are also 
2x2 matrices, corresponding to the pure state |/ifc) as- 
sociated to the measurement result of the fc-th re- 
peater, that is \r k ) — Yli,j(M rk )ij\ij)- Note that the 
computational basis i and j in the previous expressions 



are the Schmidt bases for the states \ifik) and |vfe+i) en- 
tering the repeater k. Using the fact that det(AB) = 
det(A) det(B), the previous maximization gives [26| 

N 

C N =Y[2\det(<p k )\. (5) 

fc=i 

Note that |2det(<£>fc)| = 1 if and only if \ipf.) is maximally 
entangled, which proves the announced result. 

Most of the results derived in the qubit case can be 
generalized to arbitrary dimension. Let us first consider 
the one-repeater configuration. Given a state \ip) ([I]), it 
is always possible to transform in a deterministic way 
this state into a two-qubit state of Schmidt coefficients 
(Ai, 1 — Ai) by local operations and classical communi- 
cation (LOCC). This follows from the application of ma- 
jorization theory to the study of LOCC transformations 
between entangled states [23| ■ Note that the SCP for the 
two states is the same, p ok = min(2(l — Ai), 1). 

In the case of arbitrary N, an exponential decay for 
the qubit concurrence can be shown for protocols with 
one-way communication. Given an arbitrary chain, we 
consider the almost identical chain where the first state 
is replaced by a two-qubit entangled state. It is relatively 
easy to prove that the SCP decays exponentially in the 
first chain if, and only if, it does it in the second one. We 
start with the simplest one-repeater configuration. The 
quantity to be optimized reads 

Ci =2|det( ¥ >i)|supV2|det(M r¥ > 2 )|, (6) 

M V 

where, as above, ipi (^2) is the 2x2 (d x d) matrix 
corresponding to \tpi) ( | ^2) ) , while M r is a 2 x d ma- 
trix associated to the measurement outcome r at the re- 
peater. Thus, we recognize in the r.h.s. of Eq. (|6|) the 
optimal average concurrence we can obtain out of \ip 2 ) 
by measurements on one particle which correspond to 
operators of rank 2. We denote this quantity by C, by 
Ai > A2 > . . . > Xd the Schmidt coefficients correspond- 
ing to if2 and by p ok its SCP, as above. For the outcome 
r, which occurs with probability p r , ix\ > [i r 2 denote the 
Schmidt coefficients corresponding to the resulting two- 
qubit state \ip r )- With this notation, we have 

0= 2j2PrVrtiS 2 < 2V^VT^, (7) 

r 

where 

z = 5>^2 < 1 - Ai, (8) 

r 

where the last inequality follows from the majorization 
criterion (23l |. The optimal value is obtained for x = 
p ok /2, which is achieved when p\ = 1. Thus, we obtain 

C\ = 2\det(<p 1 )\^Jp° k {2-p° k ). (9) 
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Lattice Percolation Threshold Probability 
Square i 
Triangular 2 sin ( § ) » . 3473 

Honeycomb 1 - 2 sin (■§) ~ 0-6527 

TABLE I: Bond Percolation Threshold Probabilities for some 
examples of 2D lattices. 



Note that \Jp ok (2 — p ok ) < 1, with equality if and only 
if p ok — 1. Note also, and this is important for what 
follows, that the optimal strategy only depends on |<^ 2 ), 
and not on the first two-qubit state, |<pi). 

This strategy can be generalized to the case of N 
repeaters when the measurements proceed from left to 
right. We show this generalization for the case N = 2, 
the case of arbitrary N will immediately follow. Consider 
the measurement step in the second repeater. After re- 
ceiving the information about the measurement result in 
the first repeater r*i, R 2 has to measure his particles. For 
each value of r\, and since A is a qubit, A and R 2 share 
a two-qubit pure state, \<p Tl ) ■ Therefore, for each mea- 
surement result, R 2 is back at the previous one-repeater 
situation. The optimal measurement strategy in this case 
was independent of the entanglement of the first two- 



qubit state. Thus, up to local unitary transformations, 
the measurement to be applied in the second repeater is 
independent of n, and 

C2=C 1 \yJp° 3 k (2-p?)\, (10) 

where p^ k is defined as above for the state \<fk)- It is 
straightforward that this reasoning generalizes to an ar- 
bitrary number of repeaters, so 

N+l 

C N = 2| de%0l II \Jp?( 2 -P?)- ( n ) 

Therefore, the average concurrence decreases exponen- 
tially with the number of repeaters unless the connecting 
pure states have p ok = 1. A non-exponential decay of 
the SCP when p ok < 1 would contradict this result. 
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FIG. 2: Quantum repeaters in the ID chain. The upper figure 
shows the simplest one-repeater configuration, that is gener- 
alized below. The first step of the optimal strategy for the 
one-repeater configuration and qubits, where \tpi) — 1 1^?2 ) = 
a/AT|00) + consists of entanglement swapping at the 
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FIG. 3: Example of quantum network where entanglement 
percolation and CEP are not equivalent. Each node is con- 
nected by a state consisting of two copies of the same two- 
qubit state, \ip) = i^ 2 }® 2 . The nodes marked in (a) perform 
the optimal measurement for the one-repeater configuration 
on pairs of qubits belonging to different connections, as shown 
in the inset. A triangular lattice is then obtained where the 
SCP for each connection is the same as for the two-qubit state 
1 952) - The remaining nodes perform CEP on the new lattice. 



